Resonance and absorption spectra of the Schwarzschild black hole 
for massive scalar perturbations: a complex angular momentum analysis 
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We reexamine some aspects of scattering by a Schwarzschild black hole in the framework of com- 
plex angular momentum techniques. More precisely, we consider, for massive scalar perturbations, 
the high-energy behavior of the resonance spectrum and of the absorption cross section by empha- 
sizing analytically the role of the mass. This is achieved (i) by deriving asymptotic expansions for 
the Regge poles of the S'-matrix and then for the associated weakly damped quasinormal frequencies 
and (ii) by taking into account the analytic structure of the greybody factors which allows us to 
extract by resummation the physical information encoded in the absorption cross section. 
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I. INTRODUCTION 

Since the pioneering paper of Matzner |l|inspired by a 
related but unpublished work of Hildreth |2| , wave scat- 
tering and absorption by black holes is a topic which 
has been extensively studied due to (i) its mathematical 
interest because scattering theory is a branch of mathe- 
matical physics (see, e.g., Ref. [3|) which has found with 
black holes and curved spacetimes a new and rich field of 
activities, but also and above all due to (ii) its physical 
interest in connection with various fundamental or ex- 
perimental aspects of classical and quantum gravity such 
as perturbation theory of black holes, gravitational wave 
theory, quasinormal modes and resonant scattering the- 
ory, weak and strong lensing, superradiance and instabil- 
ities, Hawking radiation, information paradox, hologra- 
phy and CFT correspondence, higher-dimensional field 
theories, analogue models of gravity... We refer to the 
monograph of Futterman, Handler and Matzner 0] and 
to references therein for the literature on wave scatter- 
ing and absorption by black holes prior to 1988 and to 
Refs. [fHHI for a short and non-exhaustive list of papers 
on this subject published since this date and which, to 
our opinion, shed light on it from new and interesting 
points of view. 

Of course, in the context of scattering and absorption 
by black holes, physicists have been mainly concerned 
with massless field theories with spin 0, 1/2, 1, 2 which 
are considered to be much more relevant, from a phys- 
ical point of view, than massive ones. However, regu- 
larly since the seventies, some interesting articles dealing 
with massive field theories and considering, in particular, 
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the influence of the mass parameter on various aspects 
of scattering in the frequency domain (scattering reso- 
nances, bound states, cross sections, instabilities...) or 
in the time domain (late-time tails, instabilities...) have 
been published (see, e.g., Refs. |15l - r26| ) and during the 
last decade an increasing number of papers dea ling with 
these same topics appeared (see, e.g., Refs. fi~2l |27H44| ). 

In this paper, we shall consider the simple case of a 
massive scalar field propagating on the Schwarzschild 
black hole and we shall revisit two particular aspects of 
scattering by this gravitational background in the frame- 
work of complex angular momentum (CAM) techniques, 
i.e., in other words, by using the Regge pole machinery. 
(For a review of CAM techniques in scattering theory, 
we refer more particularly to the monograph of Newton 
[j| and for the use of these techniques in the context of 
black hole physics, we refer to Refs. 0, 0, Q 1451450}). 
More precisely, we intend to describe the high-energy be- 
havior of the scattering resonance spectrum and of the 
absorption cross section by emphasizing the role of the 
mass parameter. It should be noted that these two top- 
ics have been already considered by numerous authors 
PUSH |H, SI, H| but precise descriptions have 
been only obtained from purely numerical analysis (see, 
however, Ref. [33| where an analytically solvable toy- 
model has been considered in order to understand the 
mass-dependence of the scattering resonance spectrum 
and a very recent paper by Hod [44| where the fun- 
damental resonances of near-extremal Kerr black holes 
due to massive scalar perturbations are derived analyti- 
cally) . The CAM approach will permit us to go beyond 
numerical considerations and to understand analytically 
in term of the mass parameter well-known effects such 
as the migration of the complex quasinormal frequen- 
cies and the behavior of the absorption cross section. Of 
course, it is important to recall that Regge pole tech- 
niques are formally valid for "high" frequencies. As a 
consequence, they do not allow us to describe in the usual 
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CAM framework the existence of the bound state spec- 
trum 0, [2(| HH HH, El| which is an important aspect 
of the massive scalar field theory on the Schwarzschild 
black hole appearing for rather low frequencies. 

Our article is organized as follows. In Sec. II, we shall 
first construct high-frequency asymptotic expansions for 
the Regge poles associated with the massive scalar field. 
We shall use two different approaches, both emphasizing 
the role played by the black hole photon sphere: (i) a 
powerful approach developed recently by Dolan and Ot- 
tewill in Ref. [47[ which is based on a novel ansatz for 
the Regge poles and the associated Regge modes and (ii) 
a more traditional approach (see Refs. [48|, |49| for previ- 
ous applications to massless theories) based on the WKB 
method developed a long time ago by Schutz, Will and 
Iyer (see Refs. (5ll - [55| ) to study the resonant behavior 
of black holes. Then, from the Regge poles, we shall ob- 
tain the weakly damped quasinormal frequencies of the 
massive scalar field in term of the mass parameter. This 
will allow us to show explicitly that when the mass of 
the scalar field increases, the oscillation frequency of a 
quasinormal mode increases while its damping decreases. 
In Sec. Ill, we shall consider the absorption problem for 
the massive scalar field. From Regge pole techniques, we 
shall make a resummation of the absorption cross sec- 
tion and then provide a simple formula describing very 
precisely, at high energies, its behavior and emphasizing 
more particularly the role of the mass parameter and of 
the black hole photon sphere. We shall finally conclude 
this paper by briefly considering some possible general- 
izations of the present work. Throughout this paper, we 
shall use units such that h = c = G = 1 and assume 
a harmonic time dependence exp(— icut) for the massive 
scalar field. 



II. REGGE POLES FOR THE MASSIVE 
SCALAR FIELD AND ASSOCIATED COMPLEX 
QUASINORMAL FREQUENCIES 

A. Generalities and notations 

In this subsection, we shall describe the first problem 
we intend to solve by using the CAM machinery and we 
shall also fix the main notations used in our article. 

We first recall that the exterior of the Schwarzschild 
black hole of mass M is defined by the metric 

ds 2 = -(1 - 2M/r)dt 2 + (1 - 2M/r)- 1 dr 2 + r 2 da\ (I) 



where da 2 



de 2 



2 denotes the metric on the 



unit 2-sphere S 2 and with the Schwarzschild coordinates 
(t,r,9,tp) which satisfy t €] — oo,+oo[, r g]2M, +oo[, 
6 e [0,?r] and <p € [0,2tt]. 

In order to simplify discussions below and to inter- 
pret physically some of our results, it is also necessary 
to point out various aspects of scattering of massive (and 
massless) particles by the Schwarzschild black hole linked 
more or less directly with the existence of its photon 



sphere at r = 3M (see, e.g., Chap. 25 of Ref. [56| or, 
for more precisions, Chap. 3 of Ref. [57( and, for some of 
the notations we shall use, Ref. [12] )• We consider a par- 
ticle with rest mass /i and energy uj > \x and we denotes 
by p(oo) — 'ui 2 — fx 2 and by v(uj) — p(u)/u) the particle 
momentum and the particle speed at large distances from 
the black hole. We note, in particular, that 
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(2) 



We recall that, associated with the parameter lu and for 
/j, fixed, there exists a sphere located at r = r c (ui) with 



r c (uj) = 2M 



1 



8v 2 {u)) 



1/2 > 



f + (l + 8v 2 (lu)) 1/2 



(3) 



on which the massive particle can orbit the black hole on 
unstable circular (timelike) geodesies. We have r c (uj) g 
]3M, 4M[. We also recall that the critical radius r c (uj) 
defines a critical impact parameter 



b c (ui) 



M 



V2v 2 (uj) 



iv 4 {ui) + 20v 2 (lj) - 1 



+ (1 + 8v 2 (lj)) 



3/2 



1/2 



(4) 



The black hole captures any particle sent toward it with 
an impact parameter b < b c (u>) while particles with im- 
pact parameter b > b c {u) are scattered. As a conse- 
quence, for particles with rest mass fj, and energy lu, the 
geometrical cross section of the Schwarzschild black hole 
IS (Tg CO (u)) = 7r6^(w) or reads more explicitly 



Cgeo(w) 



ttM 



f2 



2v i {u) 



8v 4 {uj) + 20v 2 (lo) - 1 



+ (1 + 8v 2 (lo)) 



3/2 



(5) 



It should be noted that v(ui) = I for fi = and from 
Eqs. ©-([5]) we recover the existence of the so-called black 
hole photon sphere located at r c (u>) = 3M (the place on 
which the massless particles can orbit the black hole on 
unstable circular null geodesies). Moreover, in this case, 
the corresponding critical impact parameter is given by 
b c (oj) — 3y/3M and, as a consequence, the geometrical 



cross section a 



goo 



21-kM 2 of this black hole for 



massless particles is also recovered. It should be also 
noted that for /i ^ we have the asymptotic expansions 



r c (w) = 3M 



f 



JJ_ 

9w 2 



(A 



b c (oj) = 3V3M 



1 



oVcoM = 27ttM 2 



v 

+ 3w 2 



° — 



(6a) 
(6b) 
.(6c) 



So, at high energy, a massive particle behaves as a mass- 
less one and, in particular, orbits the black hole on un- 
stable circular geodesies very near the photon sphere. 
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From now on, we consider a massive scalar field $ with 
mass /i propagating on the exterior of the Schwarzschild 
black hole. It satisfies the wave equation (□ — ^ 2 )<I> = 
which reduces, after separation of variables and the 
introduction of the radial partial wave functions ^/(r) 
with uj > and I € N, to the Regge- Wheeler equation 



drl 



+ [uj 2 ~Vi(r)] <^ = 0. 



(7) 



In Eq. (J7]), Vg(r) denotes the Regge- Wheeler potential 
given by 



V t (r) = 1 - 



2M 



a (l+l/2) 2 -l/4 2M 

M + -5 + 



(8) 

while r* is the so-called tortoise coordinate defined from 
the radial Schwarzschild coordinate r by dr/dr* = (1 — 
2M/r). Here, it is important to recall that the function 
r* = r*(r) provides a bijection from ]2M, +oo[ to ] — 

CO, +00 [. 

In this paper, we shall focus on the Z/V-modes (see, e.g., 
Ref. [58| or Chap. 30 of Ref. [5^|) which are the solutions 
of with a purely ingoing behavior at the event horizon 
r = 2M, i.e., which satisfy 



M r ) ~ T t {uj)e~ l " r * 

— > — oo 

and which furthermore, at spatial infinity r 
an asymptotic behavior of the form 



(9a) 
-oo, have 



p{u) 



1/2 



x ^[p(u)r, + (Mn 2 /p(u>)) ln(r/M)] 

+i?f(w)e +i[p(w)r * +(MAl2/p( " ))ln( '' /M)]N ) .(9b) 



In Eq. (|9bl) . p(w) which now denotes the "wave number" 
is given by 



p(u) = (uj 2 - /i 2 ) 



2\l/2 



(10) 



while T(.(bS) and Rt(u>) are transmission and reflection 
coefficients linked by 



\Ri(u;)\ 2 + |T,( W )| 2 



1 Vw > and VI e N. 



(11) 



It is interesting to note that this relation can be derived 
from the properties of the Wronskian 



,i{r),<t>u,t>{r)] = (f>u,,e( r ) ( J^T^M 



■^rMr) ) Mr)- (12) 



Indeed, from ([7]) we can show that this Wronskian is 
a constant and by evaluating it for r* — > — oo and for 
r* — > +oo taking into account the boundary conditions 



(|9]), we then obtain relation (ITT1) . It should be also 
noted that the coefficient [uj/p^)] 1 ^ 2 in Eq. (|9b[) has 
been introduced in order to simplify the form of this re- 
lation. It is moreover important to remark that the IN- 
modes are naturally and without any ambiguity defined 
by the boundary conditions ([§]) for uj > /i. However, for 
< uj < fx, the situation is a little bit more complicated 
(see Chap. 30 of Ref. [59j for more precisions): indeed, it 
is hrst necessary to go into the complex u) plane and to 
carefully take into account the various branch cuts asso- 
ciated with the functions p(io) (the two cuts ] — oo, 
and ] + (i, +oo[ along the real u) axis) and w 1 / 2 (e.g., a 
cut emanating from the origin and along the negative 
imaginary uj axis) allowing us to deal with these multi- 
valued functions; then, by working on the Riemann sheet 
in which Imp(w) > (the first Riemann sheet in the fol- 
lowing), we can define p(ui) for < uj < [i and we have in 
particular p(uj) which is a pure positive imaginary; and, 
finally, the boundary conditions Q can now be used even 
for < uj < [i. 

Moreover, it is worth recalling that the transmission 
coefficients Tg(uj) permit us to construct the greybody 
factors (the absorption probabilities by the Schwarzschild 
black hole for scalar particles with energy uj and angular 
momentum £). They are given by 



= \T e (uj)\ 2 



(13) 



and, for the massive scalar field considered here, the black 
hole absorption cross section can be expressed in terms 
of them in the form 



Cabs(w) 



•]T(2£+1)I>H. 



(14) 



Furthermore, even if we do not intend to analyze all the 
aspects of scattering by this black hole and to consider 
in particular its partial elastic cross sections, its scatter- 
ing amplitude, its differential cross section, etc. it is 
interesting to recall that all these concepts are built from 
the S'-matrix which is defined by its diagonal elements 



St(u) = (~lY +1 R e (uj). 



(15) 



Finally, it should be noted that the IN-modes alone do 
not permit us to construct a basis of the solutions of 
the wave equation but they are sufficient in order to 
describe the resonance and absorption spectra of the 
Schwarzschild black hole. In other words, here it is not 
necessary to consider the usual t/P-modes or to work 
with the alternative OUT- and DOWTV-modes (see, e.g., 
Chap. 4 of Ref. ^ or Chap. 30 of Ref. 

Let us now consider the Regge- Wheeler potential Vi (r) 
given by ([5]). From now on, we shall assume that the 
condition 



2M/i G]0,l/2[ 



(16) 



is satisfied. We shall therefore restrict our study to scalar 
field mass y, rather weak. This assumption allows us 
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to simplify our study because it saves us from compli- 
cated discussions on the behavior of Vi(r) for different 
values of the angular momentum I. Indeed, it is well- 
known (2(| [43[ that under (|16p . Vg(r) always presents 
three extrema W € N. Thanks to Tartaglia and Car- 
dano, it is furthermore easy to prove that these extrema 
rncgM, r max {£) and r min (£) satisfy r ncg {£) < < 8A//3 < 

' max 

(£) < 4M < r m i n (£) and are given by 



^neg(^) = 
^*max (^) 
^min(^) ' 



1(1 + 1) 

3M(j? 

1(1+1) 
3M/x 2 

£(£ + 1) 
3M/i 2 



2\ ^cos 



V{£) 



where 



V(£) 



3M 2 H Z 



2W-^cos 



[£ 2 {£ + lf 



m 2n 
3 3 

3 3 



(17a) 
,(17b) 
(17c) 



-9M 2 fi 2 £{£ + 1) + 9MV] (18a) 



and 



arccos 



27 



(18b) 



with 



m = 27W ^ + " 2 ™ V ^ + 1)2 

+27M 2 n 2 £{£ + 1) - 216MV] . (18c) 

[For i = 0. (fT7|l agrees with Eq. (13) of Ref. Of 
course, only the extrema r max (€) and r min (^) which lie in 
the physical region r > 2M govern the behavior of the 
Z/V-modes defined above by the Regge- Wheeler equation 
([7]) and the boundary conditions (J9j) and, furthermore, 
it is very important to keep in mind that r max (£) cor- 
responds to the peak of a local potential barrier while 
^min(^) denotes the minimum of a local potential well. 
Finally, it is interesting to note that 



WW = 3M 



1 



1 - 27M 2 /i 



2, ,2 



9(£+l/2) 2 



O 

£+l/2-> 



1 



(^ + l/2) z 



(19) 



and therefore, for large angular momenta, the peak of the 
Regge- Wheeler potential Ve(r) lies very near the photon 
sphere of the Schwarzschild black hole located at r = 3M. 

The process which has allowed us to define the IN- 
modes for < ui < /i moreover permits us to extend for 
complex ui values the diagonal matrices T, R and S. We 
recall that a pole of T^{oj) [let us note that it is also a 
pole of Re(uj) and Si(uj)] such that Ti(lu) / Ri(lu) remains 
regular and which lies in the lower half plane of the first 



Riemann sheet associated with the multivalued function 
p(u>) is a resonance of the scalar field. We also recall that 
resonances are symmetrically distributed with respect to 
the imaginary uj axis. We shall focus our attention on 
those lying in the fourth quadrant of the considered Rie- 
mann sheet. It is well-known that they can be sepa- 
rated into two families (see, e.g.. Chap. 30 of Ref. [59( or 
Refs. [16|, [20j, I23, [25|, |26|, |36|, |43| ) corresponding respec- 
tively to complex frequencies uj satisfying Re a; < /j and 
Rew > /i and therefore respectively associated with [see 
Eq. ©]: 

- a bound state spectrum (these modes are nor- 
malizable, purely ingoing at the horizon and fall off ex- 
ponentially at spatial infinity), 

- a quasinormal mode spectrum (these modes are 
not normalizable, purely ingoing at the horizon and 
purely outgoing at spatial infinity and they oscillate at 
both boundaries). 

From a physical point of view, the existence of the 
bound state spectrum is directly related to the presence 
of the potential well close to r min (l) (see also Refs. pHHS 
HI, [H, and Chap. 30 of Ref. [5j| for more precisions) 
while the existence of the weakly damped quasinormal 
modes is due to the presence of the potential barrier at 
r miH (/) or, equivalently, can be semiclassically described 
in terms of "surface waves" lying close to the photon 
sphere at r = 3M. In the remaining of this section, we 
intend to discuss more precisely this last point by using 
the CAM machinery. We are not be able to provide an 
analogous description for the bound state spectrum. 

To conclude this subsection, we shall briefly recall some 
aspects of CAM techniques we shall extensively use in 
the following. We first note that, for ui > 0, the ma- 
trices T, R and S previously defined can be analyti- 
cally extended into the CAM plane: we transform the 
ordinary angular momentum £ into a complex number 
A = £ + 1/2 and we construct the analytic extensions 
T\-i/ 2 (u), R\-i/2(u) and S , A -i/ 2 (w) of T t {uj), Rt(u) 
and Se(uj) from ©, ©, © and (15]). Let us also re- 
call that the Regge poles are defined as the poles of the 
T-matrix [or, equivalently, as the poles of the matrices 
R or S] for which 7\_ 1 / 2 (w)/-RA-i/2( a; ) remains regu- 
lar and that they lye in the first and third quadrants of 
the CAM plane symmetrically distributed with respect 
to its origin. The Regge poles can be also considered as 
"eigenvalues" associated with the so-called Regge modes 
^u,A-i/2( r ) which satisfy ([7]) and are purely ingoing at 
the horizon and purely outgoing at spatial infinity [see 
Eq. (j9])]. In Ref. [48] we showed that, for a massless 
scalar field propagating on the Schwarzschild black hole, 
the complex frequencies of the weakly damped quasinor- 
mal modes can be obtain analytically from the so-called 
Regge trajectories, i.e., from the curves traced out in the 
CAM plane by the Regge poles as a function of the fre- 
quency uj. Mutatis mutandis, the reasonin g le ading to 
these results (see for more precisions Refs. [43, |48j and 
Appendix A of Ref. [4!|) can be repeated in the con- 
text of a massive scalar field theory. Let us denote by 
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uj ln = u j£ - iT in /2 with t e N and n e N\ {0} the com- 
plex quasinormal frequencies lying in the lower half plane 
of the first Riemann sheet associated with the multival- 
ued function p(to) and by A„ with n € N \ {0} the Regge 
poles lying in the first quadrant of the CAM plane. If we 
describe the associated Regge trajectories by the func- 
tions A n = X n (uj) with lu > fj,, we have the semiclassical 
relations 



Re 



A„(^)=^+l/2 



and 



Tjn _ Im X n (uj) 
2 



d/dw ReA n (w) 



(20a) 



(20b) 



B. Regge poles: The Dolan-Ottewill method 

In this subsection, we shall obtain high-frequency 
asymptotic expansions for the Regge poles by using and 
extending to the massive scalar field a new and powerful 
method introduced and developed by Dolan and Ottewill 
in Ref. ]47|. It is important to note (see below) that, in 
this method, the critical parameters (|3|) and ((4]) associ- 
ated with the massive scalar particle are the main ingre- 
dients permitting us to construct the Regge modes and 
the Regge poles associated with the field theory. 

Extending the reasoning of Dolan and Ottewill to the 
massive scalar field, we introduce the following ansatz to 
describe the Regge modes (f> U j.\ n (u)-i/2( r ) and the corre- 
sponding Regge poles X n (uj): 



4>u,\ n (u)-i/2(r) = ?VA„( w )-i/2(?")exp 



iu>v(uj) 



2Mb c {uj) 2 /r c (uj) 2 



2 \ 2 



r c (u>) 



dr' 



(21a) 



with 



i=i j=i 



r c (w) 



exp ^Tf'^r)^)]-* 
\fc=o / 



(21b) 



and 



> , \ >fnl/ / \ \ r / m *(n), , \ \ X\ n \uiv(uj)) X\' l> (ljv(uj)) X\ rl> (u>v(uj)) X\ n> (lov(u>)) 



(n), 



>(<»)] 



MP 



>(<*)]* 



We invite the reader to compare our Eqs. (|2Tal) . (|21b|) and with Eqs. (5), (39) and (38) of Ref. [47]. It should 
be noted more particularly that, in order to extend the Dolan-Ottewill method for the massive scalar field, we must 
consider as the natural asymptotic parameter the "momentum" u)v{uS) — p(uj) instead of the energy u) and assume 
that u) > fx. 

By inserting Eqs. (I21a|) . (|21b[) and (|2"2"|) into ([7]) with ((5J) where £ —> X n (ui) — 1/2 and after a tedious calculation, we 
obtain 



X { i\ujv(uj) 



X2 (ujv(lo) 
X^ (ujv(lo) 
X^^v^) 



= 3%/3M 

= ia(n) 
_ 3V3 

.a(n) 
= *~M 2 ~ 
_ 3\/3 

~ ~w 

,a(n) 
= % ~M 4 



1 

1 + - 



3 \[wu(w)] 2 / 27 \[wu(a;)] 4 y (i/[u»(u)]->o 



1 



(9 



9 V.[^w(w)] 2 / 18 V[wu(w)] 4 / u/[^(w)]->o 
60a(n) 2 - 29 -372a(n) 2 + 115 



;(«)]< 



11664 



104976 



wv(w)] 2 / /i/[ui)(w)]->o \ [wd(w)] 4 



-1220a(n) 2 + 1357 4444a(n) 2 - 3707 



419904 



1259712 



O 



-2357520a(n) 4 + 4630008a(n) 2 - 99373 



O 



29386561536 n/[uv( u )]^o \ [ujv{uj)] 2 

144920784a(n) 4 - 439855800a(n) 2 + 28395953 



2115832430592 



O 



(23a) 
(23b) 
(23c) 
,(23d) 
(23e) 
(23f) 
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and finally 



X n (uj) — 3voM um(uj) + ia(n) + 



5 , s2 29- 3888MV 

— a(n) — 

36 v ' 432 



r i 

V (3V3M wu(w)) 



+«a(n) 



305 
3888 



a(n) 2 



1357 - 46656M V 



15552 



1 



49115 , 192917 - 4339008MV 



-a(n) 



839808 v ' 1679616 
99373 - 32192640MV + 2176782336MV 



(3y/3M uv(u)) 2 
a(n) 2 

1 



40310784 



(3V3Muw(w)) 3 



+ia(n) 



3019183 . , 4 18327325 - 311008896MV , , 2 

< n ) Tonnooo.o «( n ) 



60466176 v ' 120932352 
28395953- 6226336512MV + 117546246144MV 



2902376448 



(3V3A/ww(w)) 4 



1 



Mu»(u)->+oo \(3\/3M w;;(w)) 5 , 
In the previous equations, we have 

a(n) = n-l/2 for neN\{0}. 



(24) 



(25) 



We can also convert the previous asymptotic expansion in uv(ui) into an asymptotic expansion in lu. From (0) or 
([TO]), we have 



A„(w) = 3v3Mw + ia(n) + 



5 , , 2 29 + 1944MV 
36 aW 432 



1 



(Sv^Mw) 



-Ma(n) 



305 



3888 



a(n) 2 + 



1357 - 46656MV 



15552 



1 



+ 



49115 



192917 - 1189728M V , 2 99373 + 4339008MV + 952342272M> 



2, ,2 



r4,,4 



839808 
+za(n) 



a(n) 4 + 



1679616 



-a(n) — 



40310784 



(3V3Mw) 3 



3019183 , , 4 18327325 - 54867456A/V 



60466176 v ' 120932352 

28395953 + 611380224MV - 11 7546246 144M VI ( 



2902376448 



O 



V (3\/3Mu>) 4 J m^+oo \ (3\/3Mu;) 5 



(26) 



r 



For pi = 0, formula (|26[) is in agreement with Eq. (38) 
and Eqs. (40)-(45) of Ref. 0. 



C. Regge poles: The WKB approach 

In this subsection, we shall derive again the high- 
frequency asymptotic expansion (I26|) for the Regge poles 
but, now, we shall use a more traditional approach (see 
Refs. [ill HH for previous applications to massless theo- 
ries) based on the WKB method developed a long time 
ago by Schutz, Will and Iy er !5lTl54| (and extended to 
higher orders by Konoplya j55j) to study the resonant 
behavior of black holes and to determine more particu- 
larly their weakly damped quasinormal frequencies. We 
shall thus check formula (|26p but we shall also see that 



the Dolan-Ottewill method is, in the context of the Regge 
pole determination, a much more powerful approach than 
the WKB one providing more quickly the same results 
when we need to capture higher-order terms in asymp- 
totic expansions. Indeed, in order to derive ([25)1 we shall 
now start with a fifth-order WKB approximation for the 
greybody factors (fT3l) and, therefore, we shall work with 
very heavy expressions (see below) . Of course, if we only 
need the leading order or the next-to-leading order of the 
asymptotic expansion (|25)> . it seems to us that the WKB 
approach remains more tractable. 

Before to begin the technical part of our work, it is also 
interesting to note that the Dolan-Ottewill method devel- 
oped in the previous subsection and the WKB approach 
we shall use here are based on different physical concepts. 
Indeed, as already previously noted, the Dolan-Ottewill 
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ansatz is constructed from the critical parameters (J3)) and 
© associated with the massive particle while the WKB 
calculation will extensively use the maximum (fT7b|) of the 
Regge- Wheeler potential defining the scalar field theory. 
Of course, in both approaches, for very high frequencies, 
it is the photon sphere at t c (uj) — 3Af and the corre- 
sponding impact parameter b c (ui) = 3^/3M which play 
the crucial role. 

For I £ N and oj > with w 2 near the peak Vz(r maiX (£)) 
of the Regge- Wheeler potential, we can use, following 



Iyer, Will and Guinn [s3 - l53 | and taking into account 
some results displayed in Ref. [55|, a fifth-order WKB 
approximation for the greybody factors (|T3l) . After a 
tedious calculation, we obtain 

W = 1 + expUn] (27) 

with 



<w x ,1/2 fl 2 ^15 l2 3 \ 4 /I155 j4 315, 2l 35,, 35 , 2 5, N 



{2 " \b4 " lb / \2U48 " 2t>b b4 V2H 62 J 

/255255, 6 225225, 4 , 15015,,, 45045, 2l2 3465,,, 3465, , , 315, , 

- 63 6364 H 6365 H 6364 bib 6 636465 H 6367 

V 131072 3 32768 3 4096 3 8192 3 4 2048 3 1024 512 

1155,, 315,, 315, 2 35 , \ 8 / 66927861, H 20369349, fi , 2909907,., 

b\ + 6 4 6 6 + br b 8 h n + b% 6?6 4 + b%b b 

2048 4 512 4 6 1024 5 256 V V 8388608 3 524288 3 131072 3 5 

14549535, 4 , 2 765765, 4 , 765765,,, , 45045,,, 765765, ,,, 135135, 2 , , 

H bibj b% b 6 bib 4 b 5 H bib 7 bibi H bib 4 b 6 

262144 3 4 65536 3 16384 3 J 8192 3 7 32768 3 4 8192 3 b 

135135,2,2 9009,2, 135135, , 2 , 9009, , , 9009, , , 693 , , 45045, 4 

H 636^ H 636465 636467 6 3 6 5 6 6 H 6369 H 67 

16384 3 J 4096 3 8192 4 2048 2048 1024 32768 4 

9009, 2 , 9009, , 2 693 , , 693 , , 693 , 2 63 , \ in l 

6:6 6 646= H 6 4 6 8 -I 6567 -I bi 610 Zn > 

4096 4 4096 5 1024 8 1024 5 7 2048 6 512 / J 

,-1/2 lY 7. 3 \ / 1365,4 525,2, 95,, 85 , 2 25 \ 2 

+7rfc 1/2 <^ 63 H 64 + 6i H 6364 6365 6: H 6 6 

W 64 3 16 7 V 2048 256 3 64 3 5 128 4 32 7 

/ 285285, fi 315315,4, 28875,,, 79695,,,, 9765, 2 , 8505, , , 1365, , 

+ b% H 6^64 bib b b\b\ H 6?6 6 H 6 3 6 4 6 5 6367 

V 655 36 3 16384 3 4 2048 3 5 4096 3 4 1024 3 6 512 3 4 5 256 3 7 

2625,, 1155, , 1085,2 24 5, \ 4 / 121246125, 8 43648605, fi , 7912905,., 

H 63 6 4 6 6 bi -\ 6 8 z„ 4 + 6S H 6S64 6S65 

1024 4 256 512 5 128 / \ 4194304 3 262144 3 65536 3 

37011975,4,2 2777775,4, 2477475,,, , 225225,3, 2327325, 2 , 3 585585, 2 , , 

6,62 H 6f6 6 -I 6^6465 6^67 H bibi 6^6 4 6 6 

131072 3 4 32768 3 8192 3 5 4096 3 16384 3 4 4096 3 

555555,2,2 63525,2, 525525, , 2 , 54285, , , 49665, , , 7035, , 165165 

0,0c; -I o,Os 636/65 H 636467 H b-ibc.be, 636a 

8192 3 5 2048 3 8 4096 3 4 5 1024 3 4 7 1024 3 5 6 512 3 9 16384 

47355,2, 45045, , 2 5985, , 5355, , 5145, 2 945, \ R l 

H 6:6 6 H 6 4 6. 6 4 6 8 6 5 6 7 bi H 610 Zn ? 

2048 4 2048 J 512 8 512 5 7 1024 6 256 J °J 

,_ 3/3 f / 119119 153153 l4l 16107, 3l 47229 l2l2 6405 l2l 6237 l , , 1155, , 

+nk 3/2 < b% 6f6 4 H b%b 5 -\ 6^64 bib 6 636465 H 6367 

\V 131072 3 32768 3 4096 3 8192 3 4 2048 3 1024 512 



A 



V 131072 3 32768 3 4096 3 8192 3 4 2048 3 1024 512 

1995,, 1095, , 1107,2 315, \ /156165009, 8 63864801, fi , 13216203,. 

bl + 6 4 6 6 H bi 6 8 + 6§ 6§6 4 H 6^ 

2048 4 512 b 1024 5 256 / \ 8388608 3 524288 3 131072 3 

62777715,4,2 5354349, 4 , 4945941,,, , 519057,3, 4669665, 2 , 3 1399167 

I 636: 6 4 6 6 636465 H 6367 6 2 K H f 

262144 3 4 65536 3 16384 3 8192 3 32768 3 4 8192 

1368675,2,2 185661,2, 1279047, , 2 , 161133, , , 154917, , , 28077, , 

I 6, Or o,o« H 630,165 63O4O7 63050ft -I 63OQ 

16384 3 5 4096 3 8192 3 4 5 2048 2048 1024 3 9 

400785,4 14 3241, 2 , 139293, , 2 23457, , 22029, , 21777, 2 5607, \ 2 
•-^TeS- 54 " ^096- 5A " ^096- 54&5 + T024 6468 + WM^ 7 + ~ ^iT 610 ) ^ 



8192 3 b 



.(28) 



Here, we use the notations and 

« = = ^j^W, (29) 6 ^ W = ^|^t fOT P>2 ' (31) 

k = k[i) = -\vP{l), (30) 
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with 



V (£) ee Vi(r.)| r . = 



,{(■)] 



V«(r)U 



.(«) 



(32) 



and 



rf" 



for p > 2. (33) 



.{<)] 



It should be noted that, in Eq. ([2"5]). the term /c 1 / 2 .? 2 cor- 
responds to the first-order WKB approximation. Adding 
the terms /c 1 / 2 ^ and fc _1//2 permits us to constructed 
the second-order WKB approximation. Adding further- 
more the terms fc 1 / 2 ^ and k~ x ^ 2 z 2 permits us to obtain 
the third-order WKB approximation. Finally, adding 
the terms /c 1 / 2 2g, fc -1 / 2 Zg and fc~ 3 / 2 provides us with 
the fourth-order WKB approximation. This fourth-order 
WKB approximation can be found in Ref. [54| (see 
Eq. (12) of that paper) and is confirmed by our own cal- 
culations, even if some coefficients in Eqs. (6) and (7) 
of Ref. [54fl are wrong and if some misprints are present 
in Eq. (4). In Eq. (|2"5|) we have also taken into account 
the terms fc 1 / 2 ^ , fe -1 ' 2 ;s§ and k~ 3 ^ 2 z 2 corresponding 
to the fifth-order WKB approximation. They have been 
obtained from Ref. (5f| . It should be noted that we have 
been able to construct also the sixth-order WKB approx- 
imation of the greybody factors (f2"Tl) by taking into ac- 



count the terms 



3/2 4 
Zn 



and k 5 / 2 . 



Because we do not need it here, we do not display its 
very long expression but we can provide it upon request. 

Even if formula (|2"8")l has been obtained for £ £ N and 
u> > 0, it can be used in the complex frequency plane or, 
as we intend to do here, in the complex angular momen- 
tum plane. We consider that uj > but we transform the 
angular momentum £ appearing in the previous equations 
into the complex variable A = £ + 1/2 and we then con- 
sider the analytic extension Tx-i/2(^) of the greybody 
factors Tg(ui) defined by (|13l) and (l27|) as well as the an- 
alytic extension 5^-1/2 ( w ) °f the "phase" Sg(u}) given by 
We have 



r A -i/2H = — 



exp[25 A _i/ 2 (w)] 



(34) 



and the (Regge) poles of the greybody factors are the 
solutions A ra (w) of the equation 

5 A _i /2 (w) =i(n- 1/2)tt with neN\{0} (35) 

(here we only consider those of the poles lying in the 
first quadrant of the CAM plane) . By inserting (I17b[) into 
(|2"51) - ([3"3")1 and considering the transformation £ — >• A— 1/2, 
we obtain from (1351) 



{ujM) z = — A 



27 



2a(n) 
27 



A + 



29 - 276a(ra) i 



5832 



+i 
+ 



-1357a(n) + 1220a(n) 3 \ 2a(n) , r2 , 



209952 
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8545 - 736824a(n) 2 - 1193520a(n) 
544195584 

/6835519a(n) - 20410824a(n) 3 - 25861584a(n) E 



1 

A 

-7-6Qa(n) 



486 



M 2 /i 2 + M> 



4 ,,4 



39182082048 



- O — 

A^+oo \ A 4 



-125a(n) - 380a(n) 
7776 



M 2 fi 2 + 3a(n)AfV 



(36) 



We can then solve (l36l) perturbatively to recover (|26l . 

D. Regge poles: "Exact" versus asymptotic results 

It is now interesting to check the accuracy of the high- 
frequency asymptotic expansions previously obtained. In 
Ref. [45[, in order to determine numerically, for mass- 
less field theories, the Regge poles of the Schwarzschild 
black hole, we adapted the powerful method developed by 
Leaver 61] to calculate the complex quasinormal frequen- 
cies of the Schwarzschild and Kerr black holes (actually, 



we implemented a slightly modified version of it due to 
Majumdar and Panchapakesan and based on the Hill de- 
terminant 62]). But, here, we are confronted with a mas- 
sive scalar field and Leaver's method must be modified. 
In Ref. [36j], Konoplya and Zhidenko have shown how 
to achieve the corresponding modifications and, in this 
paper, we have combined the Leaver-Konoplya-Zhidenko 
method with the Hill determinant approach to deal nu- 
merically with the Regge pole spectrum. 

Figure Q] exhibits the exact Regge trajectories numer- 
ically calculated and permits us to show clearly the role 
of the mass parameter \x. We also observe numerically 
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Mass dependence for the first Regge pole 
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Mass dependence for the first three Regge poles 

















/l=0.4(2M=l) 
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^=0.0(2M=1) 







Mass dependence for the first three Regge poles 



li=0A (2M=1) 
(1=0.2 (2M=1) 
pi=0.0 (2M=1) 



2 4 6 8 10 

Re A„ (<a) 

FIG. 1: Regge trajectories for the massless and the massive 
scalar fields (small masses satisfying US}). The Regge poles 
A„(cj) (n = 1,2,3) are followed for 2A/cj = 0.03 -¥ 4. For 
more readability, the plot of Ke\„(uj) is displayed for the 
only first Regge pole. 



that the Regge poles are only defined for lj > fi. 

We have also compared the exact Regge poles numeri- 
cally obtained with the results provided by the series (|26|) 
or by the series (|24p . The accuracy is very impressive for 
high frequencies (i.e. for |2Mw| > 2) and the errors for 
both series are of the same order. However, it should be 



TABLE I: A sample of Regge pole values for the scalar field: 
Exact versus asymptotic results for (i = (2M = 1). The 
series (|24p is always truncated at the smallest term. 







Exact 


Exact 


Asymptotic 


Asymptotic 








Tin \ ( / , 1 1 




Tin A i / 1 1 1 
1111 y\ n yiu j 


1 


0.5 


1 282821 


W p 'J -L 'J O U t: 


1 284SQ8 


KJ • <J -L <J \J <J \J 










( n 19 °7«\ 






1.0 


2.586845 


0.504736 


2.586891 


0.504737 










(-0.0018 %) 


(-0.0002 %) 




2.0 


5.190075 


0.501236 


5.190077 


0.501236 










(-0.000039 %) 


(0. %) 


2 


0.5 


1.452818 


1.430687 


1.469628 


1.379357 










(-1.2 %) 


(3.6 %) 




1.0 


2.688355 


1.479587 


2.690232 


1.477585 










(-0.07 %) 


(0.14 %) 




2.0 


5.242994 


1.494926 


5.243088 


1.494880 










(-0.0018 %) 


(0.0031 %) 



TABLE II: A sample of Regge pole values for the scalar field: 
Exact versus asymptotic results for (i — 0.2 (2M = 1). The 
series (124[) is always truncated at the smallest term. 







Exact 


Exact 


Asymptotic 


Asymptotic 


n 




Re A„(cj) 


Im X„(oj) 


Re A„(w) 


ImA„(u) 


1 


0.5 


1.245875 


0.506823 


1.250055 


0.504331 










(-0.34 %) 


(0.49 %) 




1.0 


2.569229 


0.502503 


2.569344 


0.502477 










(-0.004 %) 


(0.005 %) 




2.0 


5.181378 


0.500679 


5.181381 


0.500679 










(-0.000066 %) 


(0. %) 


2 


0.5 


1.411681 


1.405188 


1.4155982 


1.345471 










(-0.3 %) 


(4.2 %) 




1.0 


2.670066 


1.473100 


2.671532 


1.471040 










(-0.055 %) 


(0.14 %) 




2.0 


5.234201 


1.493274 


5.234282 


1.493228 










(-0.0015 %) 


(0.0030 %) 



noted that they increase slightly with the order of the 
Regge pole and with the magnitude of the mass param- 
eter fi. For rather low frequencies, the series (|2~4"]) is the 
most accurate one. TablesH HlIII present a sample of Regge 
pole positions calculated from this asymptotic formula. 
They exhibit a good agreement for the first Regge pole 
even for low frequencies. For the second Regge pole (and 
also for the third one not displayed here), the accuracy 
remains correct for |2Mw| > 1. 



E. Effect of the mass on the complex quasinormal 
frequencies 



Formulas (|20a|) and (|20bl) permit us to construct an- 
alytical approximations for the resonance excitation fre- 
quencies and the damping Ti n /2 of the quasinormal 
modes. By inserting into these two semiclassical formulas 
the series (|26p truncated after the term in 1/(3V3Muj) 2 , 
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TABLE III: A sample of Regge pole values for the scalar field: 
Exact versus asymptotic results for fi — 0.4 (2M = 1). The 
series (|24[1 is always truncated at the smallest term. 







Exact 


Exact 


Asymptotic 


Asymptotic 






Rp \ d,A 
rv,e A n {Lu j 


1111 A n \ LU ) 


rv,e An \yJ j 


1111 A n [us ) 


1 


0.5 


1.114489 


462801 


1.077877 


42Q4SQ 










yo.o /o) 






1.0 


2.514597 


0.495198 


2.514479 


0.495068 










(0.005 %) 


(0.026 %) 




2.0 


5.155079 


0.498976 


5.155077 


0.498975 










(0.00004 %) 


(0.00020 %) 


2 


0.5 


1.258631 


1.286677 


1.55611 


0.983336 










(-24 %) 


(24 %) 




1.0 


2.613157 


1.451939 


2.611183 


1.450803 










(0.08 %) 


(0.08 %) 




2.0 


5.207607 


1.488218 


5.207596 


1.488185 










(0.00021 %) 


(0.00219 %) 



we obtain 



,(o) _ 



3V3M 



(£+1/2) 



5 . , 2 29 + 1944MV 

a(n) 2 H — 

36 w 432 



1 



(£+1/2) 



(9 



1 



£^+00 V + 1/2) 3 
r<?n _ a{n) ( 
2 ~ 3^3A/ I 



(37a) 



235 , . 2 313- 116640MV 



3888 



O 



a(nY 



15552 



1 



1 



\ (£ + l/2Y 



(£+l/2) 2 

(37b) 



It should be noted that, for /i = 0, formula sJ37p are in 
agreement with the results obtained in Ref. [53| (see also 
Ref. [3). 



From Eq. (|37|) . we can note that when the mass // of 
the scalar field increases, the oscillation frequency of a 
quasinormal mode also increases while its damping de- 
creases. This well-known result observed numerically by 
various authors (see, e.g., Refs. [26j, |36J) is here analyt- 
ically described. It is moreover possible to provide a 
dual explanation of this result by considering the physi- 
cal interpretation of the Regge poles in terms of "surface 
waves" located close to the photon sphere 0, EE EE El] , 
even if we think that this description must taken with 
a grain of salt. In this context, we recall that ReA„ (ui) 
represents the azimuthal propagation constant of the n- 
th surface wave while Im A„ (lo) is its damping constant. 
From ([2TJ)) . it is then obvious that, for a fixed value of 
the frequency w, all the azimuthal propagation constants 
and all the damping constants decrease when the mass 
parameter fj, increases. As a consequence, the resonance 
excitation frequencies u>^ which are the frequencies for 
which a constructive interference due to the surface waves 
occurs [see also Eq. (|20ap ] must necessarily increase with 
the mass parameter /x. Furthermore, because the atten- 
uation of the n-th surface wave decreases when the mass 
parameter fi increases, the energy it radiates away dur- 
ing its repeated circumnavigations around the black hole 
also decreases; it is then natural [see also Eq. (|20b[) ] to 
observe a similar behavior for the damping of the associ- 
ated quasinormal modes. 

We can greatly improve the asymptotic expansions 
(13"T1) by going beyond the semiclassical approach. Indeed, 
the complex quasinormal frequencies are solutions of the 
equation 

A„(w) =1+1/2 with neN\{0} and I G N. (38) 

Such a result can be understood if we recall that a factor 
in 1/ cos[7rA„(a;)] appears in all the residues series over 
the Regge poles constructed from the CAM machinery 
(see, e.g., Eq. (8) of Ref. Q). Equation can be 
solved perturbatively and we then obtain 



J 



,(0) 



1 



and 



839808 



O 



3V3M 
a{n) 4 
1 



(£ + 1/2) 



36 



a{n) 2 



29 + 1944M> 



2,,2T 



432 



(£+1/2) 



17795 , sA . 18763 - 14346720MV 



1679616 



a{n) 



82283 + 20015424MV - 136048896MV 



r4,,4n 



40310784 



(£+l/2) 5 



Ijn _ a(n) f 
2 ~ 3\/3M \ 



235 , , 2 313 - 116640MV 
-a{n) 2 + — 



3888 



234857 



60466176 



a(n) A — 



15552 



653125 - 953881920AfV 
120932352 



{1+1/2)' 



a(nf 



{I- 



4832407 - 3269372544MV + 29386561536MV 
2902376448 



1 



(£ + 1/2)' 



1 



<-++«, V(^+ l/2) e 



1/2)3; 
(39a) 



(39b) 
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For /i = 0, formula (l39l is in agreement with Eq. (7) and 
Eqs. (17)-(22) of Ref. E3- 



III. CAM DESCRIPTION OF THE 
HIGH-ENERGY ABSORPTION CROSS SECTION 
AND ROLE OF THE MASS 

A. Generalities 

In this section, we shall focus our attention on the ab- 
sorption cross section defined by (fTl| (see also Refs. [H, 
[43j for previous numerical investigations). We shall more 
particularly provide a simple and very accurate approxi- 
mation of this series emphasizing the role of the mass pa- 
rameter /1 and permitting us to describe qualitatively and 
quantitatively its behavior for high frequencies/energies. 

With this aim in mind, it is necessary to replace the 
sum over the partial waves (|14|) by the series over Regge 
poles 



Cabs(^) = O-gco(w) 



bH] 2 sinKA„M-l/2)] j 

+ ... (40) 

In Eq. (|40[) . er gco (w) is the geometrical cross section of the 
black hole given by ([5]) and the A ra (w) with n £ N \ {0} 
are those of the (Regge) poles of the analytic extension 
T\_i/2(u>) of the greybody factor Te(u>) lying in the first 
quadrant of the complex A plane while the 7 n (w) are the 
associated residues. We can obtain (|4U|) from (TT4"|) by 
using the CAM machinery and by repeating for to > u, 
mutatis mutandis, the main steps of Sec. II of Ref. [14fl 
where we considered the absorption cross section for a 



massless scalar field. This is possible due to the prop- 
erties of the Regge poles and of the analytic extension 
of the T- matrix mentioned in Sec. II. A. In Eq. (|40[) . the 
dots are associated with a background integral along the 
imaginary A axis. We shall assume that it can be ne- 
glected numerically and physically at high energies. 

It is now possible to consider the high-energy behavior 
of the absorption cross section (fT4|) by replacing into (j40|) 
the functions A ra (u;) and 7 n (w) by their high-frequency 
asymptotic expansions. We have already at our disposal 
the expansion (|26| for the Regge poles A„(w). It re- 
mains to us to construct the analogous expansions for 
the residues 7 n (w). 



B. Regge-pole residues of the greybody factors 

It should be first noted that the Dolan-Ottewill method 
we have found very efficient for the construction of the 
high-frequency asymptotic expansion of the Regge poles 
cannot be naturally adapted to the determination of the 
residue asymptotic expansions (see also a remark below) . 
However, it is possible to achieve such a job very effi- 
ciently from the WKB approach of Sec. II. C. 

The Regge poles are the poles of the analytic extension 
F\-i/2{u) of the greybody factors given by (|34| and the 
solutions of Eq. ([35]) . As a consequence, it is easy to 
prove that the corresponding residues are given by 



7«M 



■1/2 



[d«S A _i/a(w)/dA]| A=sAn((i>) 



(41) 



Then, from (gTJ, (23} , (TT7b|) and <j26j), we obtain the 
asymptotic expansion (a fifth-order WKB approxima- 
tion) 



7„(w) = -—+ia(n) 

Z7T 



36tt 



305 2 . . 1357 - 46656MV 
-a (n) 



2592tt 



31104tt 



+ia(n) 



49115 



419904tt 



a 2 (n) 



15095915 



192917 - 1189728MV 
1679616tt 



18327325 - 54867456MV 



4 ^ uu .^v» u g 2 

1209323527T v ; 80621568?r v ' 

28395953 + 611380224A/V - 11 7546246 144MV 4 
58047528967T 



(3V^A/a 



O 



tUJ) 



(42) 



It is interesting to note that, in a very recent article tunately, this requires lot of work. 
[o3j . Dolan and Ottewill have developed new techniques 
based on their ansatz in order to construct analytically 
the black hole excitation factors associated with the com- 
plex quasinormal frequencies. Their techniques could be 
adapted to obtain the two first terms of (|42|) but, unfor- 
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C. CAM approximation for the absorption cross 
section 

Let us now insert (|26|) and |42|) into (|40[) and take into 
account Eq. (|Bc|) . By noting that the contribution of the 



Regge poles with n > 1 is practically negligible and that, 
for mathematical coherence, we only need the first three 
terms of (|26|) and the first two terms of (|42l) . we obtain 



<r abs (w) ps 27ttM 2 1 



v 

3w 2 



87re~ 



[27r(3\/3Af)a;] 



167re 



_ 27r sin [4tt(3\/3M)c 



2 7 r(3v / 3M)w 4tt(3v / 3M)u; 
4 7 r 2 e - 7r [-39 + 7tt + 9727rAf V] cos [27r(3\/3Ai>] \ 



27 



[2tt(3a/3M) 



(43) 



r 



Exact and CAM descriptions of the absorption cross section 




11=0.30 (2M=1) CAM approximation 

— ^=0.15 (2M=1) CAM approximation 

(1=0.30 (2M=1) Exact 

(1=0.15 (2M=1) Exact 



0^ 
(1.0 



1.0 



FIG. 2: The exact absorption cross section (|14[1 and its CAM 
approximation (|45|l are compared for /i = 0.15 and p = 0.30 
(2M = 1). The limiting constant value 27ttM 2 corresponding 
to the geometrical cross section of the photon sphere is also 
indicated. 



For n = 0, formula fl33) reduces to Eq. (2.12) of 
Ref. [64| and is the superposition of an eikonal contri- 
bution (see also Ref. [lj]) and a fine structure. It is 
interesting to recall that the eikonal contribution is the 
sum of the geometrical cross section of the black hole 
photon sphere and a sine function involving the geomet- 
rical characteristics (orbital period and Lyapunov expo- 
nent) of the null unstable geodesies lying on this photon 
sphere. It describes accurately the high-energy behav- 
ior of the absorption cross section and, in particular, of 
its regular and attenuated oscillations around a limit- 
ing value (see, e.g., Ref. [65j for a numerical analysis of 
this behavior obtained a long time ago). The fine struc- 
ture is a slightly more complicated function which also 
involves the geometrical characteristics of the null un- 
stable geodesies lying on the photon sphere and which, 
above all, permits us to capture small fluctuations lying 
beyond the eikonal contribution. For (i ^ 0, even if it is 



not so natural, an analogous description can be provided: 
we can associated the first three terms of (|4"3")l. i.e., 



= 27nM' I 1 



Sire 



_^sin [2tt(3V3M) 
2tt(3V3M) 



(44) 



with the eikonal description and the two last ones, i.e., 



27ttM 



27I sin [47r(3\/3M) 



+ 



47r(3V3Af) 

A^e- 71 [-39 + 7tt + 9727rM 2 /i 2 ] cos [27r(3\/3Af)o 

[2tt(3V3M) 



27 



(45) 



with the fine structure. We can observe that the mass 
parameter fx appears in both contributions, greatly mod- 
ifies the eikonal part and slightly corrects the amplitude 
of one of the terms (the smallest one for very high fre- 
quencies) of the fine structure. 

We have tested numerically formula (|4"3"|) . It describes 
very accurately the absorption cross section (jT4")) as well 
as its oscillations around the geometrical cross section 
and, in particular, it takes into account very correctly 
the contributions of the mass parameter fx. In Fig. [51 for 
two values of the mass parameter [i, we have displayed 
the exact absorption cross section numerically obtained 
from (TTJI by solving the problem defined by ([7]), flSJ and 
© and we have compared it with the result provided 
by (|43[) . The agreement is remarkable even for rather 
low energies. We can however observe that, for low ener- 
gies, the accuracy of the CAM description decreases when 
the mass parameter fi increases. Moreover, it should be 
noted that our numerical absorption cross sections seem 
to be in agreement with those of Ref. [43[ but totally 
disagree with those displayed in Ref. [Uj]. 

In Fig. |31 for the same values of the mass parameter //, 
we have also tested the accuracy of the CAM description 
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Exact and CAM descriptions of the fine structure (2Mi(=0.15 and 2M=1) 




CAM approximation 

Exact 



Exact and CAM descriptions of the fine structure (2Mi(=0.30 and 2M=1) 




FIG. 3: Fine structure for fi = 0.15 and /j, = 0.30 (2M = 
1). We compare the exact fine structure (|46l) with its CAM 
approximation (|45l) . 



(|45|) of the exact fine structure defined by 

_fluct,fi: 

\W) = cr abs 
( 2/i 2 



Anuci. nnc / \ / \ 

°abs ' M = ^absM 



27,Af-|l + ^-8^ Sin[2 - ( y M ^ 
y 3cj 2 27r(3\/3M)w 

(46) 

and numerically calculated. The agreement is remarkable 
for high frequencies and remains robust even for rather 



low ones. It should be noted that the behavior of the fine 
structure is very simple for 2Mfi = 0.15 because, in this 
case, the coefficient of the second term in Eq. (|45|) can 
be neglected. In fact, it vanishes for 2Mfi s» 0.149265... 



IV. CONCLUSION AND PERSPECTIVES 

In the present paper, by working in the CAM plane, 
we have emphasized explicitly the role of the mass pa- 
rameter in the resonance and absorption spectra of the 
Schwarzschild black hole for massive scalar perturbations 
and simplified considerably their description. Our work 
could be extended to the case of the massive scalar field 
propagating on the Reissner-Nordstrom, Kerr and Kerr- 
Newmann black holes as well as to the more interesting 
case, from a physical point of view, of massive fermions 
19]. This last problem could even have nice applica- 
tions in the context of multimessenger high-energy as- 
trophysics. With this aim in view, the recent articles by 
Dolan, Doran, Lasenby and coworkers [IH, |37], H|| would 
constitute a natural and solid starting point. 

It is finally important to point out that, in this pa- 
per, the Regge pole machinery has not permitted us to 
understand the existence of the bound state spectrum 
E IM HI HE HI [M IH associated with the massive 
scalar field theory defined on the Schwarzschild black 
hole. This is really a pity and we consider that it is 
one of the crucial problems to be tackled in the frame- 
work of CAM techniques. Indeed, if we could solve such 
a problem, we could maybe provide a dual simple expla- 
nation of the realization by massive fields of the so-called 
"black hole bomb" scenario (6(| or, in other words, of the 
instabilities induced by the bound state spectrum in the 
presence of rotation and due to the sup erradiance phe- 
nomenon (see, e.g., Refs. 0, H HH, M, HUE HI El). 



Acknowledgments 

We are grateful to Alexei Starobinsky for information 
concerning the Russian literature and for having pointed 
to our attention Refs. [2(| and (23|. We also thank Sam 
Dolan for providing us with recent references on black 
hole instabilities. 



[1] R. A. Matzner, J. Math. Phys. 9, 163 (1968). 

[2] W. W. Hildreth, Ph. D. Thesis (Princeton University, 

Unpublished, 1964). 
[3] R. G. Newton, Scattering Theory of Waves and Particles 

( Springer- Verlag, New York, 1982), 2nd ed. 
[4] J. A. H. Futterman, F. A. Handler, and R. A. Matzner, 

Scattering from black holes (Cambridge University Press, 

Cambridge, 1988). 
[5] P. Anninos, C. DeWitt-Morette, R. A. Matzner, 



P. Yioutas, and T. R. Zhang, Phys. Rev. D 46, 4477 
(1992). 

[6] N. Andersson and K.-E. Thylwe, Class. Quantum Grav. 

11, 2991 (1994). 
[7] N. Andersson, Class. Quantum Grav. 11, 3003 (1994). 
[8] S. R. Das, G. W. Gibbons, and S. D. Mathur, Phys. Rev. 

Lett. 78, 417 (1997). 
[9] J. Maldacena and A. Strominger, Phys. Rev. D 56, 4975 

(1997). 



14 



K. Glampedakis and N. Andersson, Class. Quantum 
Grav. 18, 1939 (2001). 

V. Cardoso, O. J. C. Dias, J. P. S. Lemos, and S. Yoshida, 
Phys. Rev. D 70, 044039 (2004). 

S. Dolan, C. Doran, and A. Lasenby, Phys. Rev. D 74, 
064005 (2006). 

L. C. B. Crispino, S. Dolan, and E. S. Oliveira, Phys. 
Rev. Lett. 102, 231103 (2009). 

Y. Decanini, G. Esposito-Farese, and A. Folacci, Phys. 
Rev. D 83, 044032 (2011). 

P. A. Collins, R. Delbourgo, and R. M. Williams, J. Phys. 
A: Math. Nucl. Gen. 6, 161 (1973). 

N. Deruelle and R. Ruffini, Phys. Lett. B 52, 437 (1974). 
N. Deruelle and R. Ruffini, Phys. Lett. B 57, 248 (1975). 
T. Damour, N. Deruelle, and R. Ruffini, Lett. Nuovo 
Cimento 15, 257 (1976). 

W. G. Unruh, Phys. Rev. D 14, 3251 (1976). 

I. M. Ternov, V. P. Halilov, G. A. Chizhov, and A. B. 

Gaina, Izv. Vuz. USSR Fiz. 21, 109 (1978). 

T. M. Zouros and D. M. Eardley, Ann. Phys. 118, 139 

(1979). 

S. Detweiler, Phys. Rev. D 22, 2323 (1980). 

L. A. Kofman, Phys. Lett. A 87, 281 (1982). 

A. B. Gaina, Sov. Phys. JETP 69, 13 (1989). 

O. B. Zaslavskii, Class. Quantum Grav. 7, 589 (1990). 

L. E. Simone and C. M. Will, Class. Quantum Grav. 9, 

963 (1992). 

H. Koyama and A. Tomimatsu, Phys. Rev. D 63, 064032 
(2001). 

H. Koyama and A. Tomimatsu, Phys. Rev. D 64, 044014 

(2001) . 

R. A. Konoplya, Phys. Lett. B 550, 117 (2002). 

L. Xue, B. Wang, and R.-K. Su, Phys. Rev. D 66, 024032 

(2002) . 

L. M. Burko and G. Khanna, Phys. Rev. D 70, 044018 
(2004). 

E. Jung and D. K. Park, Class. Quantum Grav. 21, 3717 

(2004) . 

A. Ohashi and M. Sakagami, Class. Quantum Grav. 21, 
3973 (2004). 

H. Furuhashi and Y. Nambu, Prog. Theor. Phys. 112, 
983 (2004). 

M. J. Strafuss and G. Khanna, Phys. Rev. D 71, 024034 

(2005) . 

R. A. Konoplya and A. V. Zhidenko, Phys. Lett. B 609, 
377 (2005). 

C. Doran, A. Lasenby, S. Dolan, and I. Hinder, Phys. 
Rev. D 71, 124020 (2005). 

A. Lasenby, C. Doran, J. Pritchard, A. Caceres, and 



S. Dolan, Phys. Rev. D 72, 105014 (2005). 
[39] V. Cardoso and S. Yoshida, J. High Energy Phys. 07, 
009 (2005). 

[40] R. A. Konoplya, Phys. Rev. D 73, 024009 (2006). 

[41] R. A. Konoplya, A. V. Zhidenko, and C. Molina, Phys. 

Rev. D 75, 084004 (2007). 
[42] S. Dolan, Phys. Rev. D 76, 084001 (2007). 
[43] J. Grain and A. Barrau, Eur. Phys. J. C 53, 641 (2008). 
[44] S. Hod, Phys. Rev. D 84, 044046 (2011). 
[45] Y. Decanini, A. Folacci, and B. P. Jensen, Phys. Rev. D 

67, 124017 (2003). 
[46] Y. Decanini and A. Folacci, Phys. Rev. D 79, 044021 

(2009) . 

[47] S. Dolan and A. C. Ottcwill, Class. Quantum Grav. 26, 
225003 (2009). 

[48] Y. Decanini and A. Folacci, Phys. Rev. D 81, 024031 

(2010) . 

[49] Y. Decanini, A. Folacci, and B. Raffaelli, Phys. Rev. D 

81, 104039 (2010). 
[50] S. Dolan, L. A. Oliveira, and L. C. B. Crispino (2011), 

arXiv:1105.1795 [gr-qc]. 
[51] B. F. Schutz and C. M. Will, Astrophys. J. 291, L33 

(1985). 

[52] S. Iyer and C. M. Will, Phys. Rev. D 35, 3621 (1987). 
[53] S. Iyer, Phys. Rev. D 35, 3632 (1987). 
[54] C. M. Will and J. W. Guinn, Phys. Rev. A 37, 3674 
(1988). 

[55] R. A. Konoplya, Phys. Rev. D 68, 024018 (2003). 

[56] C. W. Misner, K. S. Thorne, and J. A. Wheeler, Grav- 
itation (W. H. Freeman and Company, San Francisco, 
1973). 

[57] S. Chandrasekhar, The Mathematical Theory of Black 

Holes (Oxford University Press, Oxford, 1983). 
[58] D. G. Boulware, Phys. Rev. D 11, 1404 (1975). 
[59] B. S. DeWitt, The Global Approach to Quantum Field 

Theory (Oxford University Press, Oxford, 2003). 
[60] V. P. Frolov and I. D. Novikov, Black Hole Physics 

(Kluwer Academic Publishers, Dordrecht, 1998). 
[61] E. W. Leaver, Proc. R. Soc. London A 402, 285 (1985). 
[62] B. Majumdar and N. Panchapakesan, Phys. Rev. D 40, 

2568 (1989). 

[63] S. Dolan and A. C. Ottewill (2011), arXiv:1106.4318 [gr- 
qc]. 

[64] Y. Decanini, A. Folacci, and B. Raffaelli, Class. Quantum 

Grav. 28, 175021 (2011). 
[65] N. Sanchez, Phys. Rev. D 18, 1030 (1978). 
[66] W. H. Press and S. A. Teukolsky, Nature 238, 211 (1972). 



